Improved variational approach for the Cornell potential 
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The approximate radial wave functions for the Cornell potential describing quark-antiquark in- 
teraction are constructed in the framework of a variational method. The optimal values of the 
variational parameters are fixed by the fulfillment of the requirements of the virial theorem and 
the minimality condition for integral discrepancy. The results of calculation with the simple trial 
function are in a good agreement with exact numerical results. 
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Many investigations have been devoted to the quantum mechanical non-relativistic description of quark-antiquark 
bound states with the Cornell potential [l[ by means of different approximation procedures (see, e.g., @, H, ELUf and 
references therein) based on perturbation theory, the variational method and their combinations. Recently a new 
improved variational approach Q was proposed. This method was applied to the radial Schrodinger equation with 
^v^j , nonsingular power-law potentials r s with positive s Q. In the present work, we apply our approach to the Cornell 
£S| ■ potential with the Coulomb singularity. 

Following [l[ we study the radial Schrodinger equation in the dimensionless form 

^ ' a d 2 1(1 + 1) 
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' We use the simple functions M 
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as the trial functions, where 
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for the Cornell (Coulomb plus linear) potential 



G ltn (b, x) = x l+1 cM-^Ln^ 1 {2x b ). (4) 

21 + 1 

Here L n h (2x b ) is the Laguerre polynomial, a and b are variational parameters, N is a normalization factor (< 
ipi,n(a,b)\ipi in (a,b) >= 1). Note that these functions reproduce the exact results for the Coulomb potential if we 
choose 6=1. In order to fix the parameter values we consider two conditions to which the trial functions should 
satisfy. 

' The first is the virial theorem. The virial condition 

d 2 1(1 + 1) lk 1 
dr 2 + ~~r 2 2r ~ 2~' 



< ipi,„(a, 6)| - — + _ 2 - -- - -r|W(a, b) >= (5) 



leads to the following equation for a : 
where 

u 2 (b) 



a 3 = ku 2 (b)a 2 + u 3 (b) (6) 
JcT ' Gl n {b,x)x- 1 dx 
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The exact solution of this equation allows us to express the parameter a via the parameter b : 
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Thus, the considered problem is transformed into a one-parameter problem. Note that the virial condition is equivalent 
to the usual condition 

d < <tp l>n (a,b)\H\ijj ltn (a,b) > 

da " °- (8) 

Denote the one-parameter trial functions with a = ao(b) as 

< r|V,°,„(6) >= <Pi,„(b,r) =< r\^ n (a Q (b),b) >= ^,„(a (6), b, r). 
Then the energy has the form 

E%(b)=<*l>l n (b)\H\iPl n (b)>. (9) 

In addition to the energy, expectation values of the squared Hamiltoninan characterizing the goodness of the approx- 
imate eigenf unctions, can be calculated: 



E^{b) = ^<^ n {b)\H^l n {b)>. (10) 
As the second requirement imposed on the trial function, we select the requirement that integral discrepancy 

, m <^ n (b)\H 2 \^Jb)> 

di n ( O) = ; ; 1 , (11) 

(< ^Jb)\ H Kn(b) >) 2 

is minimum [||. The quantity di^ n (b) characterizes goodness of the approximation and is equal to zero for an exact 
solution of the Schrodinger equation. The integral discrepancy is simply connected with the local discrepancy Q 

by relation 

d l>n (b) =< D l>n (b)\D l>n (b) > . (13) 
The corresponding equation for selection of parameter b is 

^d nl (b)=0. (14) 

Note that we must find the absolute minimum of di n (b) corresponding to the minimal rotation of a trial vector under 
the action of the HamiltOnian in Hilbert space 9} . 

The minimality condition for integral discrepancy is alternative to the minimal sensitivity condition (To| 

!*e>m-a ,15) 

There can be several stationary points in the case of the function E^(b), and there is no method to select one of 
them. Besides, energy is not a preferred quantity in comparison with the other expectation values. It is also well 
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TABLE I: Linear potential (k = 0). 
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TABLE II: Cornell potential. 
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known, that calculation of the energy with increased accuracy does not always lead to the improvement of other 
characteristics. 

We compare our results with the results E^™ 1 of some variational method (VFM) [2] and with the results 
E"™ of numerical integration of the Schrodinger equation [l[ . We also compare our values of the quantity 

<v > hn =J Q [-^^) dr (16) 

with values from [l|. 

Table 1 shows that the proposed approximation gives fairly accurate results in the case of the purely linear potential. 
Table 2 demonstrates the efficiency of our approach in the case of the combined Coulomb plus linear potential. 
For example we compare our variant (14) and usual variant (15) in the case k = 1, / = 1, n = 1: 

£$(14) = 4.4619, E[ 2 l(U) = 4.4620, < v 2 > M (14) = 1.1868 (a = 0.5439, b = 1.6769), 

£$(15) = 4.4583, £$(15) = 5.0658, < v 2 > M (15) = 0.9905 (a = 1.3352, b = 0.9999) 

while the numerical solution gives E™ m = 4.4619, < v 2 > I 1 m 1 m = 1.1864 . 

We see that in spite of the simplicity of the used trial functions, the optimal choice of the variatonal parameters 
leads to the satisfactory approximation. At the same time we assume that the approximation can be improved by 
means of some complication of the trial function reproducing the behavior at the origin more accurately. 
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